Two-photon decays of highly excited states in hydrogen 
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The relativistic and nonrelativistic approaches for the calculations of the two-photon decay rates of highly 
excited states in hydrogen are compared. The dependence on the principal quantum number (n) of the ns, nd, 
and np initial states is investigated up to n = 100 for the nonresonant emissions. For the ns states together with 
CN) ■ the main E1E1 channel the contributions of higher multipoles (M1M1, E2E2, E1M2) are considered. For the 

np states the E1M1 and E1E2 channels are evaluated. Moreover, the simple analytical formula for the E1M1 
decay is derived in the nonrelativistic limit. 



I. INTRODUCTION 



The question about the role of the two-photon decays in the context of the hydrogen recombination in early Universe was 
raised at first in Hi. It was shown that the electron recombination occurs mainly into the 2s state, followed by the 2s state 
decays via two-photon E1E1 emission into the ground state. The transition frequencies of these photons are not resonant (the 
energy of each photon released will not excite any of neighbouring atom), and, therefore, the medium is transparent for these 
frequencies. The contribution of such transition rate to the microwave background was further investigated in the literature, see, 
e.g., JH-dl). The recent success in observation of the cosmic microwave background temperature and polarization anisotropy 
draws attention to the details of cosmological hydrogen recombination history. This, in turn, requires an accurate knowledge of 
q • the two-photon decay processes. 
• ' The theoretical formalism for the two-photon decay has been introduced by Goppert-Mayer ((I and the first calculation of 
the two-photon E1E1 transition 2s — > 2j(El) + Is belongs to Breit and Teller |@]. Further improvements including relativistic 
corrections were performed in works (H,[9|]. The two-photon transition 2p — Is occurring via E1E2 and El Ml channels was first 
O !* evaluated in ITolllTll . In recent paper by Amaro et al. lfl2tl these results were confirmed and further two-photon transition rates 
l— ~~ 1 ■ from initial bound states with n = 2, 3 have been evaluated. 

In connection with increased accuracy of the astrophysical experiments (see, e.g., lfl3ll ) it might become also important to 
consider atomic recombination into highly excited states in hydrogen, as proposed by Dubrovich and Grachev in J4j], and further 
(■-<■) ' developed in [ 14-16]. The magnitude of the nonresonant radiation from highly excited states is of the same order as for the 
2s — Is transition and, therefore, the same contribution to the emission escape can be anticipated. Highly excited states can 
decay by cascade emission (i.e, the transition occurs via intermediate state) and by "pure" (nonresonant) two-photon emission. 
In the case of the cascade transition the emitted photons can be absorbed immediately by the neighbouring atoms and, therefore, 
£NJ • the medium is not transparent for the resonant radiation, but it is for the nonresonant one. For this reason the question of 
separation of the "pure" and cascade emission arises. There are some papers considering this problem (see, e.g., IU7l - l20ll ). In 
the papers lfl7l[T8ll the possibility of separation between cascade and "pure" two-photon emission has been shown. However, 
in our paper lU9tl it has been demonstrated that it is impossible to divide resonant and nonresonant two-photon emission due to 
J> . the interference term between them. Recently, in [20] Jentschura accepted our point of view on the treatment of this problem. 
The separation of the cascades and the "pure" two-photon emission remains ambiguous. In works I14lllql the influence of the 
highly excited states to the Sobolev escape probability for Lyman-a photons was computed, employing the "1+1" scheme. In 



the vicinity of the Lyman-a resonance the shape of the "1+1" emission profile remains Lorentzian and the resonant/nonresonant 
or resonant/resonant interference terms do not appear. Another approach for the account of the two-photon processes from the 
highly excited states to the recombination problem was presented by Hirata in lfl5ll introducing a distinction between regions 
with resonant photon contributions and those with "pure" two-photon contributions. 

The main aim of the present paper is to evaluate rigorously the dependence of the two-photon decay rates on the principal 
quantum number of initial state n. Since the nonresonant photons being transparent for the medium contributes to the microwave 
background, it is important to investigate the n-behaviour of the "pure" two-photon emission. For this reason we calculate the 
differential transition rate for different n at the equal frequencies of the emitted photons, the energy region where the contribution 
of the cascades is almost negligible for high n. In the calculations we employ both the relativistic and nonrelativistic approaches. 
By the nonrelativistic approach we mean: (i) the long-wavelength approximation for the photon emission operator, in the case 
of dipole transitions this yields the dipole operators, therefore we will use the term "dipole approximation" for this; (ii) matrix 
elements evaluated with the nonrelativistic (Schrodinger) electron wave functions. In the relativistic approach the full photon 
emission operator together with the relativistic electron wave functions are employed. Since the orbital size of the highly excited 
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states increases quadratically with the increasing of n, the long-wavelength (dipole) approximation seems to be inadequate 
0]. The argument of the photon wave function is no longer small and, therefore, the restriction to the first term only in the 
power-series expansion of the Bessel function is under the question. One could expect that the utilization of the full photon 
emission operator contributes essential to the two-photon decay rates from the highly excited states. For the same reason the 
higher multipoles contributions (besides the dominant E1E1 channel) of the two-photon transitions for the ns and nd levels 
could become important. Moreover, since the highly excited states are strongly mixed, we should consider, the decays from the 
np states as well. 

The paper is organized as follows: In Sec. II we present the general formulas for the two-photon decay. In Sec. Ill we calculate 
E1E1 differential transition rates ns/nd — > Is at the equal photon frequencies in hydrogen atom. The influence of the higher 
multipoles (M1M1, E2E2, and E1M2) is also investigated. As next step we consider the two-photon E1E2 and E1M1 transition 
rates of the np states (Sec. IV and V). The investigation of the dependence on n is performed for all these transitions. Moreover, 
in Sec. V we derive simple analytical formula for the E1M1 decay rates of the np states in the nonrelativistic limit. 

Apart from Sec. II and III where the relativistic units are useful for convenience, the atomic units are used throughout the 
paper. 



II. GENERAL FORMULAS: TWO-PHOTON DECAY 

The two-photon transition probability A — » A' + 27 corresponds to the following second-order S'-matrix elements 

(A'\S^\A) = e 2 J d A x l d i x 2 {^ A ,{x l ) ltJ , 1 Al i {x l )S{x 1 x 2 ) ltl2 Al 2 {x 2 )^ A {x 2 )) , (1) 

where £(2:1X2) is the Feynman propagator for the atomic electron electron. In the Furry picture the eigenmode decomposition 
for this propagator reads (e.g. lulll . see also lfl9lo 

00 _ 



S(xix 2 ) = -L / i.^-^M 



- iO) + w 



where the summation in Eq. (0 extends over the entire Dirac spectrum of electron states N in the field of the nucleus, tp^ (x) is 
the electron wave function, Em is the electron energy. In Eq. 7^ are the Dirac matrices. The wave function of the photon 

A^\ X ) = J^ e M e ^--*) (3) 
V cj 

is characterized by the momentum k {lj = \k\) and polarization vector ejL (p, A = 0, 1, 2, 3), x = (r, t). For real transverse 
photons we have 

A(x) = yj^ztf^ = sJ^A^e-^ . (4) 

Integrating over time and frequency variables, taking into account photon permutation symmetry and introducing the amplitude 

U A 'A as 

S A 2 ?l = -2md{E AI +L0 + L0' - E A )XJ A 2 ?l , (5) 

yields the expression 
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where e is the electron charge, a is the vector incorporating the Dirac matrices. The labels A, A' and N abbreviate the complete 
set of the atomic electron quantum numbers (principal quantum number n, total angular momentum j, its projection m, and 
parity I) for the initial, final and intermediate states, respectively. The transition probability is defined via 



dWfj) = 2ttS{E a - E A , - u - w') 
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Setting dk = oj 2 dvdoj and integrating over oj yields 
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E N — E A + oj 



dvdv' doj' . 
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Here the summation over N represents the summation over all sets of quantum numbers of the intermediate state {nljm}. The 
sums over the projections of the total angular momentum of the final state A' and the averaging over the projections of the total 
angular momentum of the initial state A in Eq. ([8} are also included. 

Expanding further the photon wave function in a multipole series (see, e.g., Refs. 1I2III22I0 we consider contributions with 
different multipole structure E1E1, E1M1, E1E2, E1M2, Ml Ml, and E2E2. The nonrelativistic limit can be easily obtained from 
Eqs. (0]i and ®, replacing the photon and electrons wave functions by their nonrelativistic analogue. In this way we describe 
electric dipole, magnetic dipole and electric quadrupole photons (see, e.g., Ref. lf23lo . 



III. E1E1 DECAY OF THE ns, nd STATES 

First, we focus on the decay rate of the ns and nd levels (A = ns/nd — > A' = Is) in hydrogen. The common expression for 
the two-photon transition probability in relativistic evaluation can be written as 



dW A 2 2)(u) = e 
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duj , 



(9) 



where the summation over the photon polarization and the integration over the photon angles have been carried out, u' = 
E A — E A i — oj. Here Aj^ f (u;) are the multipole components of the transition operator, the symbol A indicates magnetic 
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multipoles (A = 0) or electric multipoles (A = 1). The multipole components Ay),(w) are defined in the same way as in 



21, 22, 24]. They are given in the transverse gauge by the expressions 



.(0) 



JM 
.(1) 



(oj) = jj(ujr) aY^{r) , 



A%(w) = ( J SK) + : ^)afW(r) + v / J(jTT) 



jj(ojr) 



(10) 



where jj[x) is the spherical Bessel function, Yj^(r) are the vector spherical harmonics. The transverse gauge is also called 
velocity gauge since the electric dipole matrix elements in this gauge turn into the velocity form dipole amplitudes in the 
nonrelativistic limit. In the length gauge the magnetic multipole components have the same form, whereas the electric multipole 
components are given by 
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-jj+i(ur) 3?W(t) + 



J 



J 



-j J+ i(wr) aYj A p(r) 



J+l 



J 



jj(ur)Yju(r)l, 



(11) 



where I is the identity operator, Yjm(y) are the spherical functions. Explicit formulas for the one-electron matrix elements 
Ajw(u) in the length and velocity gauges can be found in j22ll . For the case of E1E1 transition one can easily obtain 



dW 



E1E1 
AA' 



t 327TOJ oj' 
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E A 
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N 



E a +oj 



doj . 



(12) 



The nonrelativistic limit for the E1E1 transitions can be obtained by expanding the Bessel functions in the photon wave 
functions Eqs. dTOb-dTTb and keeping only the nonrelativistic value for the large component of the electron wave functions. Thus, 
we come to the following expressions in the length gauge 



dW^ s {oj) = e* 



i OJ 3 oj' 3 
27n 



\S ls . ns (oj) + S ls . ns (oj')\ doj 



(13) 
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and 



with notations 



i c . .3. ./3 

dW*lf s {w) = e 4 U |5i. lWj (o;) + S lsM {J)f du: , (14) 

1O07T 



o I \ ^ ( R 1s\r\Rnip){Rn'p\r\R ns /nd) 

<3l s , ns /nd(, w J - 2^ p p — T7^ ' ( 5) 

n'rt ^ 



{Rn>vV\Rni) = / r 3 R n ,i,{r)R nl (r)dr, (16) 
o 

where R n i(r) are the radial parts of the nonrelativistic hydrogenic wave functions, and E n i are the corresponding energies. 

The corresponding decay rate for the two-photon transitions ns/nd —> Is can be obtained by integration of Eq. dT2b [in 
nonrelativistic case Eqs. ( fOt and (TT~4b l over the entire frequency interval 

W%T = \j dWjT(u), (17) 
o 

where ujq — Ea — E' A . The gauge invariance of the transition probability was investigated in [25] for the one-photon transitions. 
The two different forms for the Efc one -photon probabilities in combination with different gauges were obtained in B26I1 . In our 
paper [11] the results of [26] were applied for the evaluation of the two-photon transition probabilities. The gauge invariance 
serves as a tool for testing the numerical evaluation procedure in both relativistic and nonrelativistic case. 

The aim of this paper is the investigation of the nonresonant two-photon decay rates from the highly excited hydrogenic states, 
since these transitions could play an important role in studies of the anisotropy of the cosmic microwave background. It was 
expected that two-photon transitions from the highly excited states could contribute essentially in the recombination process. 
Accordingly, we study the behaviour of the frequency distribution function of the transition probability of the initial state with 
the increasing principal quantum number n. Comparison between the relativistic and nonrelativistic calculations are performed 
in order to examine possible deviation from the dipole approximation. The results of numerical calculations are presented in 
TableQ]for the ns — Is two-photon E1E1 transitions. For the analysis of the n-behaviour of the "pure" two-photon decays we 
present the results in terms the normalized distribution function 



1 dW E1 ? x 

qEIEI _ 1 aVV ns,ls 3 



at the equal frequencies of the emitted photons x = lo/ojq = 0.5. The choice of the frequency is caused by the absence of 
the cascades in this region. The E1E1 two-photon transition probability W^ 1 ^ 1 = 8.229 s _1 is taken as normalization factor. 
Moreover, the normalized distribution function is multiplied by factor n? in order to emphasize the dominant behaviour of Q^l^l 
for large n values. In Table Q] we present also the corresponding quantity Q„5^g for the E1E1 two-photon nd — Is transitions. 
Again the normalization factor M^ 1 ^ 1 for the frequency distribution function is taken. 

These results reveal that for large values of n the scaling of the frequency distribution function for the E1E1 two-photon 
ns/nd — Is transitions is close to 1 /n 3 and the dipole approximation holds even for the large n. This means the decrease of the 
nonresonant two-photon transitions contribution for large values of n to the recombination dynamics of the primordial hydrogen 
in the Universe. The relativistic corrections appear to be not essential for the considering two-photon transitions from the highly 
excited levels to the ground state. The Table|I]shows also, that two-photon E1E1 ns — Is transition probabilities decrease faster 
then 1 / n 3 for the low values n and picture changes for the n s» 10 when the normalized quantity Q^l^l reaches the asymptotic 
value m 6.6. For the nd — Is ElE 1 two-photon transitions behaviour is quite different. Namely, for the n 40 the values of the 
transition rates decrease slower 1/n 3 and then, as in the case of ns — Is transitions, reach the 1/n 3 asymptotics. 

To make the picture complete the relativistic evaluations have been also performed for E2E2, Ml Ml and E1M2 transitions 
as the corrections to the E1E1 ns-level decay. Their contribution are compiled in Table HI1 and appear to be negligible. The 
Table HI1 shows that the contributions of the higher multipoles to the nonresonant two-photon emission for the ns — Is decays 
is negligible. The parametric estimation for these transition rates is (aZ) in atomic units. The values of the M1M1, E2E2, 
and E1M2 transition probabilities are in good accordance with the reported in ITEzTl . It is obvious that corresponding values of 
the transition probabilities for the two-photon emission of the nd states are also negligible and all of them can be omitted in 
astrophysical investigations. 
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TABLE I: Comparison between the results of relativistic and nonrelativistic calculations of the normalized distribution functions Q nB ,i s and 
Qnd^is at the equal frequencies of the emitted photons x — 0.5 for the interval of the principal quantum number n = [2, 100] and n = [3, 100], 
respectively. 



n 


^ElEl \ 

yns,is(,nonrel.J 


^ElEl / i \ 


^ElEl \ 


n E1E1 jwi i 


2 


10.35 


10.35 






3 


0.32 / 


o coo 
0.32s 


32. 21 


32.21 


J 


/.233 


/.233 


/i o m 


49.9/ 


8 


6.674 


6.765 


56.39 


56.46 


10 


6.591 


6.648 


57.97 


58.00 


21) 


0.393 


o.4yu 


AA 1 Q 

OU. 13 


An no 
00.08 


30 


6.576 


6.460 


60.51 


60.46 


40 


6.567 


6.450 


60.65 


60.60 


50 


6.562 


6.445 


60.71 


60.66 


60 


6.559 


6.442 


60.71 


60.70 


70 


6.558 


6.441 


60.73 


60.72 


80 


6.557 


6.439 


60.74 


60.73 


90 


6.556 


6.439 


60.74 


60.74 


100 


6.556 


6.438 


60.75 


60.75 



TABLE II: The relativistic values of the normalized distribution functions Q^^ii, Q^is an d Qnl^is at the equal frequencies of the emitted 
photons x — 0.5 for the interval of the principal quantum number n — [2, 100]. 



n 


, n E2E2 ( , \ 


Qns,i s (rel.) 


Q ns ,is (rel.) 


2 


1.168 x 10" 


ii 


3.019 x 10" 


ii 


2.581 x 10" 


10 


5 


1.907 x 10" 


n 


5.277 x 10" 


n 


2.962 x 10" 


10 


10 


2.967 x 10" 


n 


5.620 x 10" 


n 


2.886 x 10" 


10 


20 


3.265 x 10" 


n 


5.916 x 10" 


n 


2.861 x 10" 


10 


40 


3.342 x 10" 


n 


5.932 x 10" 


n 


2.853 x 10" 


10 


60 


3.356 x 10" 


n 


5.932 x 10" 


n 


2.853 x 10" 


10 


80 


3.361 x 10" 


n 


5.932 x 10" 


n 


2.852 x 10" 


10 


90 


3.363 x 10" 


n 


5.932 x 10" 


n 


2.851 x 10" 


10 


100 


3.364 x lO- 


li 


5.932 x lQ- 


li 


2.851 x 10" 
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IV. E1E2 DECAY OF THE np STATES 

In this section we evaluate the two-photon E1E2 transition rates np — > Is. Relativistic and nonrelativistic calculations for the 
2p — Is E1E2 transition were previously performed in papers IToll and Hill . We will follow the paper HUl . where two-photon 
E1E2 transition probability was evaluated in different forms and gauges. For reasons of simplicity in this section we present 
the evaluation of the transition rate within the nonrelativistic length form. The electric multipole operators can be written as 
(nonrelativistic limit): 



V Ek (w) 



+ 1 2w fc+ 2 
~ (2& + 1)!! 



r k Y k _ M . 



(19) 



where Yk -m is a spherical harmonic. Accordingly, the two-photon decay rate of the atomic state A with the emission of two 
electric photons can be written as (see, e.g., Ill ill ) 



dW 



EfeEfe' 
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(w,u/) 
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N 



E 

MM' raAVrt A i 
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y, (A'\V Ek {oj)\N){N\V Ek (w')\A) 

N 



En — Ea 



■uj' 



E 

N 



(A'\V Ek (uj')\N)(N\V Ek (oj)\A) 



■ UJ 



{A'\V Ek {uj')\N){N\V Ek (lj)\A) 
En — Ea + tJ 



(A'\V Ek (u)\N) (N\V Ek (u')\A) 

N 



En — Ea 



- UJ' 



5 (uj 



En — Ea 
oj' -E A + E A > ) duduj 1 . (20) 



6 



To perform the summation over the complete set of intermediate states we employ the explicit expression for the Coulomb Green 
Function B27I1 . With the aid of eigenmode decomposition of the Coulomb Green function the probability of the two-photon decay 
process in nonrelativistic limit takes the form lUlll 



2 2 ujW 

3 3 5 2 7T 



J 2 \h{uj') + I 2 (oj)f + u 2 \h{u) + I 2 {uj')Y 



duj . 



(21) 



where 



1 

7q 



dr x dr 2 r{r\e ri * gi (E A - u; n, r 2 ) 



o o 



(22) 



and 



I 2 ((jj) = -^= J J dr 1 dr 2 r\r\ e ri ? 9i{E A ~ (*>; n, r 2 ) ■ 
o o 



(23) 



The decomposition of radial part of the Coulomb Green function reads 



9i{v, r, r') 



AZ 



exp 



[L 2l+l (2r\ L 2l + 1 f 2± 



^ (21 + 1 + n)\(n + l+l-u)' 

71=0 



(24) 



where are the generalized Laguerre polynomials. The corresponding radial integrals can be evaluated analytically. Equa- 

tions (22) and (23) can be easily generalized for the case of np — > j(El) + j(E2) + Is transitions by the replacing the radial 
function 2p state, which is equal to R 2 \(r) — ^7g e ~' by the required one. 
Finally, integrating over frequencies lu yields (ujq — E 2p — E\ s ) 



W E1E2 

vv 2p,ls 



dw™ 2 > 



lo) = 1.98896 x 1(T 5 {aZf a.u. = 6.61197 x 1CT 6 s" 1 (Z = 1) 



(25) 



The Z-dependence of the W 2p i^ transition probability is also indicated. Compared with the relativistic result the relative 
discrepancy is about 0.1%. The parametric estimation of the transition rate W§^f s 2 is (aZ) 8 in atomic units. This is a 2 times 
less than the corresponding parametric estimate for the E1E1 transition probability, which is (aZ) 6 . However, with respect 
to the achieved accuracy of the astrophysical experiments, the E1E2 transition rates could be considered as a correction to the 
two-photon processes. 

In Tableimithe relativistic results for the npi/ 2 — > 7(-ET) + r y(E2) + (n — l)s transition rates are presented. Good agreement 
with the corresponding nonrelativistic results is found. Bad convergence properties of the finite basis set representation restrict 
our calculations of the transitions between neighbouring states to n = 5. This can be checked by considering the level of 
accuracy up to which gauge invariance is preserved in the numerical evaluations. For the nl — > 2j + Is transitions the gauge 
invariance is preserved with good accuracy up to the n = 100. In Table UTTl we also display the relativistic and nonrelativistic 
values for the E1E2 normalized distribution function 



Q 



E1E2 

"Pi/2, 



1 



dW E1E2 , 

np 1/2 ,ls 3 

ri 



2 W E1E2 ls dw 



(26) 



at the equal frequencies of the emitted photons £ = uj/luo = 0.5. One can see from the table that the behaviour of the nonresonant 
E1E2 two-photon emission with the increasing n values is the same as for E1E1 probability of the ns/rwi-states emission. And 
again the nonrelativistic dipole approximation works well. 



V. E1M1 DECAY OF THE np STATES 

A parametric estimate shows that the E1M1 transition probability is of the same order of magnitude as for E1E2, therefore, 
it should be included in our consideration. For the first time the E1M1 transition rate for the 2pi/ 2 — > Is emission process was 
evaluated in fllOllllll . Our results were later confirmed by Amaro et al. in B12I1 . 
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The expression for the E1M1 two-photon transition probability is given by 



M E MMmAm A i 



x (A'\V E1 (uj)\N)(N\V M1 (u;')\A) (A'\V m1 (lj')\N)(N\V e1 (lj)\A) 



x E N -E A +ui> ^ E n -E a +oj 

(A'\v e1 (u')\n)(n\v u1 (oj)\a) (A'\v m1 (u)\n)(n\v e1 (lj')\a) 2 



^at-Sa+w ^ E N -E A +U)' 



duj . (27) 



In the nonrelativistic limit the corresponding magnetic dipole operator reads V M1 (ui) — y^|/iow 3 / 2 (jiM M + ®iMm)' where 

Ho = a/2 is the Bohr magneton, jia/ m and sijvf M are the spherical components of the total angular momentum and the spin op- 
erator (spherical tensors of rank 1) of the electron. Since the operator for the magnetic photon includes total angular momentum 
and spin operator, coupled wave functions characterized by the set of quantum numbers A*" = {nlsjm} should be used, i.e. 

<Pnls jm = J2 C lZ sm s Rnl (0*$ (^Xsm, , (28) 
mim s 

where Xsm B (s = 1 /2) is the spin function. The magnetic potentials in Eq. (f27J do not depend on radial variables. Thus, only the 
intermediate states with nl = n A l A or nl = n A d A > will contribute to the transition probability in Eq. ( |27| i. Performing angular 
integrations and summations over all projections one arrives at the expressions 

dWf™)M =« 6 ^ (j) 12 u,a/ 3 du, (29) 

and 



5 12 w ° 

<™ ls 4/<">)^7(?) fufa-ufdu,, (30) 
o o 

with Lu — E 2pi/2 — Ei s . As the final result we obtain (see also ITlll ) 

W E1M1 ls = ^|^-a.u. = 9.6769 x lO"^ 1 (Z = 1) . (31) 



Again the Z-dependence of the ls transition probability is explicitly indicated. Comparison with the result of a fully 

relativistic calculation now reveals a deviation of about 0.1%. 

In case of arbitrary principal quantum number n due to the condition 

{R nl \Ml\R n , V ) ~ Snn'SlV (32) 

only the diagonal matrix elements remain in the sum over the intermediate states in Eq. d27b . Therefore, in the nonrelativistic 
limit no cascades will arise in the considering processes involving the magnetic dipole photons. We should note that, in the 
relativistic case the cascades will arise, however, the contribution of the cascades appears to be negligible in the case of hydrogen 
atom. 

For the E1M1 transition probability we have 



R np (r)r 3 Ri s (r)dr 



duj . (33) 



After the integration over radial variable r and the frequency ui within the interval [0, E npi/2 — Eu] we obtain 



1 \ 2n 

n — 1 



W E1M1 , = ( - — - ) (aZf (34) 

For n ^> 1 it follows W?„ 1M1 ,„ w — - — (l ] — — r—z, i.e. the same scaling law arises as for E1E1 transitions. 

F1/2 ' 135-7T \ 3n 2 J exp 4 n d 
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TABLE III: The relativistic values of the transition probabilities W^p™^ and ^npy 2 (n~i)s m units s 1 are presented in second and 

fifth columns, respectively. In the third and fourth columns the relativistic and nonrelativistic values for the distribution function Qnp^? 2j i s 
are given at the equal frequencies of the emitted photons x = 0.5. In the sixth and seventh columns the relativistic and nonrelativistic values 
for the distribution function Q np ^} 2> \ s are presented at x = 0.5. 



n 


w , 

nv-\ /o . (to— 


l)s 


Qnpf /2 ,is(rel.) 


✓-\P1P9 / i \ 
Qnpf/ 2 ,l S ( n0nreL ) 


W C -L JVA..L 

nth /o. f TO — l)s 


i\/ri / n \ 


s-\ P 1 TVT 1 / i \ 

Qnpf /2 ,is(nonrel.) 


2 


6.612 x 10 


-06 


15.31 


15.32 


9.682 x 10" 06 


9.99 


13.33 


3 


3.737 x 10 


-08 


3.31 


3.32 


1.189 x 10" 08 


12.65 


14.24 


4 


9.753 x 10 


-10 


14.26 


14.27 


1.989 x 10~ 10 


13.60 


14.51 


5 


6.221 x 10 


-11 


21.91 


21.91 


1.025 x 10 -11 


14.04 


14.63 


10 






34.64 


34.53 




14.63 


14.78 









20 






38.28 


38.45 




14.78 


14.82 








30 






38.94 


39.13 




14.81 


14.83 








40 






39.36 


39.22 




14.82 


14.83 








50 






39.33 


39.47 




14.82 


14.83 








60 






39.39 


39.52 




14.82 


14.83 








70 






39.43 


39.56 




14.82 


14.83 








80 






39.46 


39.58 




14.82 


14.83 








90 






39.47 


39.60 




14.83 


14.83 








100 






39.48 


39.61 




14.83 


14.84 









For the npx/ 2 — > "/(El) + j(Ml) + 2s emission process we obtain in the same way 



piiui 4 / 32 \ (aZ) 

Similarly, for n > 1 the scaling is W^ M1 2s — — 1 - ' v ' 



n Pl/2 ,2s ~ 135?r ^ 3n 2 yexp 8 n 3 ' 
In Table Hill the relativistic results for the np\i% — > "/(El) + j(Ml) + (n — l)s transition rates are presented. These results 
are in a good agreement with the corresponding nonrelativistic values, which can be easily obtained in terms of Eqs. 133l -d3Bt. 
Again we restrict ourselves to n = 5 for the transitions between neighbouring states [npi/ 2 and (n — l)s]. In Table Hill we also 
display the relativistic and nonrelativistic values for the E1M1 normalized distribution function 

at the equal frequencies of the emitted photons x = w/cjo = 0.5. As in previous cases the behaviour of the nonresonant E1M1 
two-photon emission is 1/n 3 . The nonrelativistic dipole approximation holds in this case as well. 



VI. CONCLUSIONS 



Relativistic and nonrelativistic calculations have been performed and compared. The numerical evaluation have been carried 
out employing the dual-kinetic-balance finite basis set method [ 28] with basis functions constructed from B-splines 12911 . For 
the nonrelativistic calculations the analytical expression for the nonrelativistic Coulomb Green function has been utilized. 

We have proved that the dependence of the transition probabilities on the principal quantum number of the initial state n out of 
the cascade regions is close to 1/n 3 . Therefore, the contribution of the highly excited states turn out to be much less significant 
for the astrophysical purposes as expected. Previously, the n-dependence of the nonresonant contribution was estimated in 
fl fl^ll . It was found that the nonresonant transition rates scale roughly linear, when increasing towards larger n. The nonresonant 
two-photon rates was defined, e.g., in B1411 . via neglecting the resonant state in the summation over the entire spectrum of the 
intermediate states. However, the term with the resonant state alone contributes both to the resonant rate and to the nonresonant 
one. Therefore, to our mind, the investigation of the n-behaviour of the nonresonant emission should be performed by analyzing 
the differential transition rate in the region, where the contribution of the cascades is negligible. 

The relativistic calculations were performed aiming for the search of the influence of effects beyond the nonrelativistic dipole 
approximation. The gauge invariance served as a check of the calculations. The nonrelativistic calculations were performed in 
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the length gauge. In principle, the gauge invariance also could be used in this case (see, for example, flUl ). but comparison 
with the corresponding relativistic values is enough for our purposes here. We have compared the relativistic and nonrelativistic 
results to understand whether it is necessary to go beyond the dipole approximation. It was expected, that for the highly excited 
states (n S> 1) the dipole approximation might be not very accurate due to the large argument of the Bessel function. However, 
we were able to show that even for the n = 100 the deviation between the nonrelativistic dipole approximation and relativistic 
theory is not significant. The relative difference for the corresponding values does not exceed 1.8%. Since we always consider the 
two-photon decays into the ground state, this can be the reason that the dipole approximation is valid even for such high values of 
n. The presence of the ground state provides a short-range cutoff for one of the radial variable (r, for instance, in Eq. (l24li). then 
the range of second radial variable r' will also be small, because of the exponential suppression factor exp[— (r + r')] coming 
from the Green function d24l i. Moreover, we have evaluated M1M1, E2E2, and E1M2 ns — Is transition rates. It was shown that 
they behave like l/n 3 also and the parametric estimation for them is (aZ) 10 in atomic units. 

On the score of the highly excited states are fully mixed we have considered also np — Is two-photon transitions. In the frame 
of the present work we have evaluated the E1E2 and El Ml two-photon transitions rates for the np — Is emission processes 
and investigated the behaviour of the nonresonant decay rates as function of the principal quantum number n of the initial 
state. For large values of n the nonresonant emission of the two-photon E1E2 and E1M1 emission is proportional to l/n 3 . The 
nonrelativistic dipole approximation is again in a good agreement with the relativistic calculations. Within the framework of the 
relativistic approach we have evaluated also the total decay rates for the np —> (n — l)s E1E2 and E1M1 transitions up to n = 5, 
and compared them with the nonrelativistic results. The restriction to n = 5 is set by the level of accuracy upto which the gauge 
invariance is achieved in the numerical evaluations. The relativistic calculations based on the finite basis set method does not 
converge good in this case and others methods should be employed. Moreover, we have derived analytically simple formulas 
for the two-photon np —> j(El) + 7 (Ml) + Is and np — > 7 (El) + j(Ml) + 2s emission processes within the nonrelativistic 
approach. It has been shown that in the presence of one magnetic dipole photon in the two-photon emission process the transition 
probability does not contain cascade contribution in the nonrelativistic limit. Account for the relativistic corrections leads to 
appearance of the cascade contributions. However, the smallness of the two-photon E1E2 and E1M1 transitions leaves them 
behind the astrophysical investigations. 

The main conclusion that we can stay is the rapid reduction for the all considered nonresonant two-photon transition rates 
with the increase of n. The dependence l/n 3 has been established for "pure" two-photon transition rates, corresponding to the 
radiation escape from the interaction with the matter. It was shown that the nonrelativistic dipole consideration of the two-photon 
transition rates from the highly excited states is valid even for large n, in case that the final state is one of the lowest state of the 
atom. 
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